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Abstract. In this note, we extend the results about the fluctuations of the 
matrix entries of regular functions of Wigner random matrices obtained in |23| 
to Wigner matrices with non-i.i.d. entries provided certain Lindeberg type 
conditions for the fourth moments are satisfied. In addition, we relax our 
conditions on the test functions and require that for some s > 3 
J M (1 + 2|fc|) 2s |/(fc)| 2 dfc < 



1. Introduction and Formulation of Main Results 

Let Xn — -j=Wn be a random Wigner real symmetric (Hermitian) matrix. In 
the real symmetric case, we assume that the entries 

(W N ) jk , l<j<k<N, 
are independent random variables such that the off-diagonal entries satisfy 
E(W N ) jk = 0, Y(W N ) jk = ct 2 , 1 < j < k < N, m 4 := sup E{W N )% < oo, (1.1) 

and the Lindeberg type condition for the fourth moments takes place, 

L N (e) -> 0, as N -> oo, Ve > 0, (1.2) 

where 

L "^ = jp E E (K^| 4l {|(^) y |>eViV})- (1-3) 
l<i<j<N 

Here and throughout the paper, E£ denotes the mathematical expectation and V£ 
the variance of a random variable £. 

In addition, we assume that the diagonal entries satisfy 

E(W N )u = 0, 1 < i < N, al := supE\(W N )u\ 2 < oo, (1.4) 

i,N 

ijv(e) -> 0, as N ->• oo, Ve > 0, where (1.5) 

l ^ = Jj E E {\( W ^\ 2l {\(W Nht \>,VN})- (i- 6 ) 
l<i<Ar 

We note that (|1.2p and (|1.5p are satisfied if 

sup E\(W N )ij\ 4+e < oo, su P E|(T^Ar) 4l | 2+e < oo. (1.7) 



S.O'R. has been supported in part by the NSF grants VIGRE DMS-0636297 and DMS-1007558. 
D.R. has been supported in part by the NSF grants VIGRE DMS-0636297, DMS-1007558, and 
DMS-0905988 . 

A.S. has been supported in part by the NSF grant DMS-1007558. 



2 S. O'ROURKE, D. RENFREW, AND A. SOSHNIKOV 

If {A=(W N )u, l<i<N, (W N ) jk , 1 < j < k < N, } are i.i.d. N(0, er 2 ) random 
variables, Xm belongs to the Gaussian Orthogonal Ensemble (GOE). 
In the Hermitian case, we assume that the entries 

Ke(W N ) jk , 3m(W N ) jk , l<j<k<N, {W N ) ih l<i<N, 

are independent random variables such that the off-diagonal entries satisfy 

Emt(W N ) jk = E3m{W N ) jk =0, 1 < j < k < N, (1.8) 

2 

Yme(W N ) jk = V3m(W N ) jk = 1 < j < k < N, m 4 := sup E\(W N ) jk \ A < oo, 

(1.9) 

and the Lindeberg type condition (11.21) for the fourth moments of the off-diagonal 
entries takes place. 

In addition, we assume that the diagonal entries satisfy 

E(W N )u = 0, 1 <i<N, a\ i^supEK^),,! 2 < oo, (1.10) 

i,N 

and the Lindeberg type condition (|1.5I) for the second moments of the diagonal 
entries takes place. 

If {^{W N )uA <i <N, me(W N ) jk , 3m(W N ) jk , 1 < j < k < N, } are i.i.d. 

N(0, ^-) random variables, belongs to the Gaussian Unitary Ensemble (GUE). 
We define the empirical distribution of the eigenvalues of X N as 



JV 

1 



1 N 



N 

i=l 

where Ai < . . . < Ajv are the (ordered) eigenvalues of X^. 

Wigner semicircle law (see e.g. [35], [5J, PQ, [2]) states that the random measure 
fix N (dx,uj) converges almost surely in distribution to the (non-random) Wigner 
semicircle distribution fi sc . The limiting distribution is absolutely continuous with 
respect to the Lebesgue measure and its density is given by 

^r(z) = - xH Wa] (x). (1.12) 



Its Stieltjes transform 



,,,( : . : l d -^M = Z ~^^\ * e C\[-2aM. (1-13) 



is the solution to 



a^gi(z) - zg a {z) + 1 = (1.14) 



that decays to as \z\ — > oo. 

This paper is devoted to the question of the fluctuations of matrix entries 
of f(Xpf) for regular test functions /. Lytova and Pastur ([18]) considered the 
GOE/GUE case and proved that 

N(f{X N ) %] -E{f{X N ) ij ))^N{Q,'^^J 1 {f)), (1.15) 
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with /3 = 1(2) in the GOE (GUE) case, 

^ 2 (/) := V(/(T7)) - \ f" 7 (/(*) - f(y)) 2 j^ V^ 2 - x*y/*° 2 ~ V 2 dxdy, 

J — 2tx J — 2cj 

(1.16) 

where r\ is distributed according to the Wigner semicircle law (|1 . 12|) . 

In |23j . Pizzo, Renfrew, and Soshnikov considered the non-Gaussian case and 
proved the following theorems. 

Theorem 1.1 (Theorem 1.3 in [23]). Let Xn — -jj^Wn be a random real symmet- 
ric Wigner matrix such that the off-diagonal entries (Wjv)ife) 1 < j < 
k < N, are i.i.d. random variables with probability distribution [i and the diagonal 
entries (Wn)u, 1 < i < N, are i.i.d. random variables with probability distribution 
Mi- 

Let f : R — > R be four times continuously differentiable on [—2(7 — 5,2a + 5] for 
some 5 > and h{x) be a C°°(R) function with compact support such that 

h(x) = lforx£ [-2a -S,2<7 + 5], 6 > 0. (1.17) 

Then the following holds, 
(i) For i — j, 

VN(f(X N ) vl -E((fh)(X N ) vl )) - ^1(W N ) U -> JV(0, »?(/)), (1.18) 

a 

in distribution as N — >■ oo. where 

vl(f) ■= 2 (c 2 (/) - a 2 (/) + ^/3 2 (/)) , (1.19) 



«(/):= E(/(r?)^) = z/M^^-^da:, (1.20) 

/?(/) := E (f^t^L^ = JL 1^ /(x)^ 2 - a 2 )^V4^^, (1.21) 
w 2 (/) defined in \1.16\) . and K4(/j) is the fourth cumulant of \i, 

K4 ( M ) = / ^(dx) - 3( / x 2 fi(dx)) 2 = E|(W^v)i 2 | 4 - 3cr 4 . 



/// is seven times continuously differentiable on [—2a — 5,2a + 5], then one can 
replace E((fh)(X N ) u ) in flJgj) by 

,-2<j 



f(x)-^\/4,a 2 -x 2 dx. (1.22) 
2ira 



(ii) For i ^= j, 



NifiX^ij -E((/fc)(*tf)tf)) - — (W0v)« iV(0,d 2 (/)) (1.23) 

(T 

m distribution as N — )■ oo, where 

d 2 (f):=uj 2 (f)-a 2 (f). (1.24) 

I// is six times continuously differentiable on [—2a — 5, 2a + 5], then one can replace 
E((fh)(X N )ij) in byO. 

(Hi) For any finite m, the normalized matrix entries 

VN(f(X N ) l3 - E((/fc)(*w)ii)) , 1 < * < 3 < m, (1.25) 
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are independent in the limit N — > oo. 

Remark 1.1. If f £ C 4 (R) and ||/||4,i < oo, where 

H/IIni := max ( / \d k f / dx k {x)\dx ] < oo, (1.26) 

0<fc<« / 

then one can replace E((fh)(Xif)ij) in U.18H1.23\) byE(f(X N ))ij. 

In the Hcrmitian case, the analogue of Theorem 11.11 was proved in Theorem 1.7 

of [53]. 

Theorem 1.2. [Theorem 1.7 in [55]/ Let Xm = TTjy^JV &e a random Hermitian 
Wigner matrix il.8^1.lU\) , such that the off-diagonal entries (WN)jk, 1 < j < k < 
N, are i.i.d. complex random variables with probability distribution /i and the di- 
agonal entries (Wn)u, 1 < i < N, are i.i.d. random variables with probability 
distribution fii. 

Let f : M. — > R be four times continuously differ entiable on [—2a — 6,2a + 6] for 
some 6 > 0, and h(x) be a C°°(R) function with compact support satisfying fi.lTj ). 
Then the following holds. 

(i) For i = j, 

s/N (f(X N ) u -E((fh)(X N ) u )) - —(W N ) tl -> N(0,v 2 2 (f)) (1.27) 

a 

in distribution as N — > oo, where 

vl(f):= u, 2 (f)-a*(f) + ^f3 2 (f), (1.28) 

w 2 (/),a(/), and /?(/) are defined in H1.16\) . \1.2(J\) , and \1.21}) . and is given 

by 

K 4 (p) :=E\(W N ) 12 \ 4 -2a\ 

If f is seven times continuously differ entiable on [— 2a — 8, 2a + 8], then one can 
replace E{(fh)(X N ) zl ) in \L~Wl\j by FT2~2}) . 

(ii) For i =/= j, 

VN(f(X N ) l3 -E((fh)(X N ) l3 )) - ^!l(W N ) l3 -)■ 7V(0,d 2 (/)), (1.29) 

(7 

m distribution as N — > oo, where N(0,d 2 (f)) stands for the complex Gaussian 
random variable with with i.i.d real and imaginary parts N(0 :2 -d 2 (f)) : and d 2 (f) 
defined in jl.24\ ). 

If f is six times continuously differentiable on [—2a — 6,2a + 8], then one can 
replace E ({fh){X N ) tj ) in H~2~9\) by 0. 

(Hi) For any finite m, the normalized matrix entries 



<N (f(X N ) l3 - E{{fh){X N ) l3 )) , 1 < i < j < m, (1.30) 

are independent in the limit N — > oo. 

Almost simultaneously with [23], Pastur and Lytova (see Theorem 3.4 in [22] ) 
extended the technique of [18) and proved the convergence in distribution for the 
normalized diagonal entries \fN (/(Xjv)ii — E (/(A0v)m)) , 1 < i < iV, when the 
real symmetric Wigner matrix Xn has i.i.d. entries up from the diagonal and, in 
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addition to the requirements of Theorem 11.11 the cumulant generating functions 
log (Ke zWl2 ) is entire. The results of [22] hold provided the test function satisfies 



(l + 2\k\y\f(k)\dk < co, 

i 

where f{k) is the Fourier transform 



/(fc) = -^= / C ifeX /l • (hoi) 



The approaches of [23. and [22 a are independent from each other. In particular, 
Pastur and Lytova prove the convergence of the characteristic function of 



'N (f(X N ) u -E(f(X N ) u )). 

In addition, in the non-i.i.d. case, Theorem 3.2 of [35] proves that 

W[VN (f(X N ) u - E(f(X N ) u ))} -> 2v 2 (f) 

provided the matrix entries (Wisr)ij are independent up from the diagonal and 
satisfy 

nW N )jk = 0, V(W N ) jk = a 2 , E(W N )% = m 3 ,E(W N )% = m 4 < oo, (1.32) 
sup E\(W N ) jk \ 6 < oo. (1.33) 

j,k,N 



In this paper, we extend Theorems 11.11 and 11.21 to the non-i.i.d. setting provided 
the matrix entries satisfy the fourth moment Lindeberg type conditions ()1 .2[) and 
(11.361) for the off-diagonal entries and the second moment Lindeberg type condition 
(|1.5[) for the diagonal entries. Moreover, we relax the smoothness condition imposed 
in [23] on the test function. 

Consider the space H s consisting of the functions <fi ■ R — > R that satisfy 

iuii2 — / n i oii.i\2s 



(l + 2\k\Y s \(j}(k)\' i dk < oo. (1.34) 

The result below is valid (see Remark 11.31) provided a test function / coincides 
on the interval [—2a — 5,2a + 5] with some function from T~i s for some s > 3, 5 > 0. 
Thus, roughly speaking, we require that / has 3 + e derivatives on [—2a— 5, 2a + 6]. 

We recall that C"(R) and C n ([—L, L}) denote the spaces of n times continuously 
differentiable functions on R and [— L, L], respectively. We define the norm on 
C n ([-L,L])as 

WfWcHl-L.L]) ■= max x G [-1,1], < I < n^j . (1.35) 

Theorem 1.3. LetX^ — ~^^N be a random real symmetric (Hermitian) Wigner 
matrix M.l)) , Jj.^| ) (respectively li.fflU. lDj) such that the Lindeberg type condition 
lll.ty) for the fourth moments of the off-diagonal entries and the Lindeberg type 
condition \1.5)) for the second moments of the diagonal entries are satisfied. Let 
f G "Hs, for some s > 3. Let m be a fixed positive integer, and for 1 < i < m, 
assume that the following two conditions hold: 
(Ax) 

A,7v(e) -> 0, as N -> oo, Ve > 0, (1.36) 
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where 

A,iv(e) = ^ E E {\( W Nh\%\(w Nhj \>eN^}) ; (1-37) 
r-3¥* 

(A 2 ) 

m A (i) := lim 1 ]T E|(W W )«| 4 (1.38) 

exists. 

Then the results (i)-(iii) of Theorem ] 1.1\ (respectively Theorem 1 1 . 2\) hold for the 
joint distribution of the matrix entries {y/W (f(Xpf)ij — E (/(^Ar)ij)) , 1 < i < j < 
m}, where n 4 {p) must be replaced in 11. 1 Sty by 

Ki{i) :~ m4(i) — 3cr 4 , 1 < i < m, (1.39) 

in the real symmetric case and by 

K 4 (i) := m 4 (i) - 2<7 4 , 1 < i < m, (1.40) 

in the Hermitian case. 

In addition, the following estimates for E(/(X/v)jj) take place. 

(iv) Let f : K ->■ K freZon^ io Cj(K), tfie space o/ seven times continuously 
differentiate functions with compact support, and supp(f) € [— L, L] for some L > 
0. Then there exists a constant Const\{h,a,a\,m^) depending on L, a, o\,m 4 , 
such that for 1 < i < N, 



(1.41) 

ie< / G C 8 (R), i/jen f/iere exists a constant Coros^c, cti, m 4 ) such that 

\E(f(X N ) u ) - /(a;)_L V4<r2 - (1.42) 



< ^ Il/ll8,l,+ , l<t<N, 



where 



:= max Qf (|a;| + 1)|^ (a;)^, < Z < . (1.43) 
exists a constant Consts(a, oi,m 4 ) such that 
\E(f(X N ) jk )\ < ^^^^ I/Im, 1 < j < k < N, (1.44) 
where ||/||6,i *s defined in U.26}) . 

Remark 1.2. If the distribution of the entries ofW^ does not depend on N, The- 
orem \T73\ vroves that \/~N {f{Xjsi)a — E (/(X/v)ii)) converges in distribution to the 
sum of two independent random variables ^p-Wa and N(0,2v 2 (f )) (in the Hermit- 
ian case, the second term is N(0,v 2 (f))), and fori ^ j, \JN (/(Xjv)ij — E(/(Xjv)y)) 



converges in distribution to the sum of two independent random variables ^p-Wij 
and N(0,d 2 (f)), where in the Hermitian case 7V(0,d 2 (/)) stands for the complex 
Gaussian random variable with with i.i.d real and imaginary parts N(0, \d 2 {f)). 
This is exactly the way Theorems \1.1\ and were formulated and proven in the 
i.i.d. case in [23l. 
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Remark 1.3. If f : R — > R coincides on [—2a— 5, 2a+5] with a function <fi £ H s , for 
some 5 > and s > 3, then Theorem \1.3\ holds for (/(-Xjv))t,j — E(//i(Xjv)),j, 1 < 
*j j < m , where h G C^°(R) is defined in \1.17^ . 

If one requires that the test function / satisfies the same smoothness assumptions 
as in [23], then the extension of the results of [23 to the non-i.i.d. setting mostly 
follows the outline of the proof in [23 . To relax the conditions of Thereoms 11.11 and 
ll.2l on the test functions, we improve the estimate on the variance of the resolvent 
entries (see Proposition ^. 2p . and employ Proposition 12.21 

We will denote throughout the paper by consti,Consti, various positive con- 
stants that may change from line to line. Occasionally, we will drop the dependence 
on N in the notations for the matrix entries. Typically, we consider in detail only 
the real symmetric case as the proofs in the Hermitian case are very similar. Some 
parts of the proofs that are almost identical to the arguments in the i.i.d. case will 
be only sketched. 

The rest of the paper is organized as follows. We prove several preliminary 
results in Section [2j including Proposition 12.21 Section [3] is devoted to the bounds 
on the mathematical expectation and variance of the resolvent entries. Theorem 
11.31 is proved in Section |4l Finally, we discuss Central Limit Theorem for quadratic 
forms in the Appendix. 

2. Preliminary Results 

We start with the following lemma. 

Lemma 2.1. Let Xn = ~^Wn be a random real symmetric Wigner matrix 
\1.J$ such that the Lindeberg condition for the fourth moments of the off- 

diagonal entries and the Lindeberg condition il.5\) for the second moments of the 
diagonal entries are satisfied. Then there exists a random real symmetric Wigner 
matrix Wn and a non-random positive sequence e^v — » as A — > oo such that 

E(W N ) jk = 0, Y(W N ) 3k =a\ l<j<k<N, (2.1) 

sup E(W N )% < oo, (2.2) 

E(W N ) U = 0, 1 < i < N, (2.3) 
su P E|(VM«| 2 < oo, (2.4) 

i,N 

sup (|G^vy , 1 < i,j < a) < e N VN, (2.5) 

F(W N ^ W N ) -)• 0, as A -> oo. (2.6) 

An equivalent result holds in the Hermitian case. 

Proof. It follows from (|1.2j) and f| 1 . 5|) that there exists a non-random positive se- 
quence ejy — > as JV - ¥ oo, such that 



1 



E 



Ne 



N 



(\(W N h\\ {WNhjl > eN ^ } )^0. (2.7) 
E {\( W "h\ 2l {\(w Nht \>e N VN } ) ->0. (2.8) 



N l<i<j<N 
1 \ "* b/'|/W..VJ2 



Ki<N 
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One can always choose ejv in such a way that it goes to zero sufficiently slow. Define 
Wn by truncating the entries of Wn at the level €nVN, i-e. 

(W N ) l3 = (W N ) tJ l {l{WNhj ^ eNVW} . (2.9) 
It follows from and dHJ} that 

P(W N + Wn) -> 0, as N oo. (2.10) 
Let us now fix i < j and consider the off-diagonal entry (Wjv)ij- We note that 

7^ liv := |E(W^)«| < E (l(^ivk|l { | (Wjv)4 .|> ejvVF} ) (2.11) 

< ^3754 E (l(^)^| 4l {|(^),|> £N v^}) - ( 2 - 12 ) 
7^ := E\W N \l — cr 2 = E (l(^)d%|(y^) y |>^ } ) (2-13) 

< ^ E (lW0«l%|(^|>«Wff>) • ( 2 - 14 ) 

Then we can constract (Wjv)ij as a mixture of the random variable (W}v)y with 
weight 1 — T ^ ,JV kt^t~ and some random variable a; ,• a? with weight T )± ,N + '; 3 ^ w 

V Ncn n VNejv ive 7v 

chosen so that 



ki^ivl < ewViV, (2.15) 
E(Wjv)»j = 0, (2.16) 
n\(W N ) ij \* = o 2 . (2.17) 
It follows from our construction and (|2.7[l that 

2 



E E E (i(^i 4l {i(^i> £N y]v } )^o. 

i<j<N N l<i<j<N 

(2-18) 

The diagonal case i — j can be treated in a similar way. We write 

7-< A w := |E(W0v)«l < E (l(Wiv)ii|l { , ( ^ )w |> ew VA?>) ( 2 -19) 
< ~ ^ E (l(^)«l%|(^v)«|>^}) • (2.20) 



'Ne N 

One then constructs (Wjv)« as a mixture of the random variable (Wn)u with weight 
1 — -t#^ and some random variable a; ; w with weight -4=*^- chosen so that 

VNen ' ' V Nen 

\ai, itN \<e N ^/N, (2.21) 

E(Wjv)« = 0. (2.22) 

Then 

E F + (W N )u) < ^rE (l(^)«| 2 l { |(^) ii |>e w ViV}) 0> (2-23) 

l<i<Ar w 

as TV — >• oo. 

It follows from ([2~TU|) . ([235]) . and ([2~23|) that (JH]) is satisfied. The equations 
(l2"Tj) and follow from ([235]) , ([2~T7| . and (|2T2"2"]) . The estimates J2J1J) and (|2~i|) 
follow from the construction. □ 



ON FLUCTUATIONS OF MATRIX ENTRIES 



9 



The proof of the next result is very similar to the proof of Lemma |2~T1 and is left 
to the reader. 

Lemma 2.2. Let Wn be a random real symmetric Wigner matrix U.l\) . I{1.4\ ), and 
let tl.36\) is satisfied for 1 < i < m, where m is some fixed positive integer. Then 
there exists a random real symmetric Wigner matrix Tn and a non-random positive 
sequence ejv — > as N — > oo such that 

{T N ) jk = (W N ) jk , m + l<j,k<N, (2.24) 

P((Tjv)ife = (Wjv)ifc, 1 < i < m, 1 < k < N) —> I, as N ^ oo, (2.25) 

E(T N )ik = 0, 1 < i < m, 1 < k < N, (2.26) 

•¥(T N ) ik = a 2 , i ^ k, 1 < i < m, 1 < k < N, sup Y(T N )u < oo, (2.27) 

\<i<m,N 

sup HT N )tk < oo. (2-28) 

N,i^k:l<i<m, l<k<N 

sup \(T N ) ik \ < e N N^ 4 . (2.29) 

l<i<m,l<k<N 

The next Proposition is essentially due to Bai and Yin (see e.g. [S], [5]). 

Proposition 2.1. Let Xn = "^^v be a random real symmetric (Hermitian) 
Wigner matrix (respectively lil.8Wl.lU\) such that the Lindeberg type 

condition for the fourth moments of the off-diagonal entries and the Lindeberg 
type condition \1.5)) for the second moments of the diagonal entries are satisfied. 
Then 

\\X N \\^2a (2.30) 

in probability as N — > oo. 

Remark 2.1. Bai and Yin (]5], [2]^ considered the i.i.d. case and proved the almost 
sure convergence. However, convergence in probability is enough for our purposes. 

Proof. Because of Lemma [2.11 it is enough to prove (|2.30l) for Wn- Moreover, we 
can modify Wn by making all diagonal entries equal to zero. Clearly this changes 
the norm of Wn at most by en- The proof uses the Method of Moments. It is 
enough to show that there exists a sequence fcjv, N > 1, such that 

—Z- -> oo, JS-^ -> 0, as N -> oo, 2.31 
log N log N 

where cn is the same as in Lemma 12.11 and for any constant z > la 



E 

N 



Tr((W N /VN) 2kN 



< oo. (2.32) 



The proof of (|2.32p in ([5 ) is combinatorial in nature and does not use the fact 
that the entries are identically distributed. By Markov inequality, it follows from 
([2"32"]l that 



E 1 

N 



(||Wjv/viV|| > z) < 



for any fixed z > 2a. Therefore, by Borel-Cantelli lemma, we have 

P(\\W N /VN\\ > z i.o.) = 0, 
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which together with the Semicircle Law implies that 

\\W N /sqrtN\\ -> 2a a.s. 

□ 

The rest of this section is devoted to the bounds on V[J^° f(x)fj,(dx, w)], where 
fi(dx, oj) is a random measure on (R, B) and B is the Borel u-algebra on R, provided 
one can control V[J_ 3m Tz^f^idx, oj)] for 3m 2 ^ 0. We follow the ideas of Propo- 
sition 1 in 26 and Proposition 3.5 in |15j . In particular, our computations below 
are close to those in [26], where /j,(dx,ui) was taken to be the empirical spectral 
distribution of a random matrix. 

Let (QjJ 7 ) be a measurable space, and (fi', IF' ,V) be a probability space such 
that f2' = R x f2, and J 7 ' is generated by B x J 7 . We denote an elementary outcome 
by u/ = e 1 x SI, and consider a random variable X{u>') = x. When it does 

not lead to ambiguity, we will denote the sub-algebra {R x D, D £ J 7 } by T . Let 
us denote by w), B £ B, oj £ SI, a regular conditional distribution for X given 
J 7 , i.e. 

For each B C R, B £ B, oj -> fj,(B, oj) is a version of V{X £ B\F). (2.33) 
For a.e. oj, B —> fi(B,oj) is a probability measure on (R, B). (2.34) 

Such regular conditional distribution for X always exists (see e.g. [11]). In partic- 
ular, if / : R — > C is such that 

E|/(X)| < oo, (2.35) 

then 

/+oo 
f(x)^(dx,oj) a.s. 
-oo 

The following proposition holds. 

Proposition 2.2. Let E\X\ < oo, s > i, and f £ W s , where T-L s is defined in 
(Tlffl ). Then 

V[J f(xMdx,oj)}=Y{E(f(X)\F)} 

/*oo P OO p~\~00 -1 

< Const s \\f\\ 2 s dye-yy 2 *- 1 dxV[ 3m- -n(dt t u)]. 

Jo J -oo Joo t x — iy 



where Const s is some absolute constant that depends only on s. 

Proof. Since s > i, it follows from (|1.34l) that / £ L 1 (R) which implies that / £ 
C*o(R), the space of continuous functions vanishing at infinity. In particular, (|2.35|) 
holds and E(f(X)\T) is well defined. Since E(e lkx \T), k £ R, is L 1 continuous 
family of bounded random variables, one can write 

1 f°° . 

E(f(X)\F) = -= f(k)E(e lkX \T)dk. (2.36) 

V 27T J -oo 

Then 

-i poo poo 

V[E(/(X)|.F)] = — / / f(k 1 )f(k. 2 )C(k l ,k. 2 )dk i dk 2 , (2.37) 

W J-oo J-oo 

where 

C(k u k 2 ) = Cov(E(e iklX \F),E(e ik * x \F)) . (2.38) 
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One can rewrite the r.h.s. of (|2.37|) as 



1 

2^ 



OO />OG 



/(fci)(l + 2|fc 1 |) s /(fe)(l + 2\k 2 \) s K(k 1 ,k 2 )dk 1 dk 2 , (2.39) 



oo J —OO 



where 

K(h, k 2 ) = C(h, fc 2 )(l + 2|fc 1 |)- s (l + 2|fc 2 |)- s . (2.40) 

Therefore, 

Y[E(f(X)\T)]<l-\\ff s \\K\\, (2.41) 
where ||iT|| denotes the operator norm of the integral operator 

/oo 
K(x,y)g(y)dy. 
-OO 

It follows from f|2.38D and (|2.40p that K is a non-negative definite operator. Since 
C(&i, fea) is a bounded continuous function on E 2 , the operator 7^ is trace class and 

/oo 

K(u,u)du. (2.42) 
-OO 

Thus, 

V[E(/pO|.F)] < — H/li; jf^ C(k,k)(l + 2\k\)- 2s dk. (2.43) 

Let us fix z = x + iy, y 7^ 0, and consider 3m as a function of A. Its Fourier 
tranform is given by ^ e -\ k v\~ ikx _ Therefore, 



-i poo poo 

V[E(3m(X -x- iyT x \T)\ = - / / e -^+\ k ^e l ^~ k ^ x C{k u k 2 )dk x dk 2 . 

^ J — oo J — oo 

(2.44) 

Formally, taking into account 

JV<*-^)*dx = 2w5(k 2 -ki), 



OO 

e" 



we obtain 

fOG pOO 



/oo poo poo 

V[E(3m(X-x-iy)- 1 \J r )]dx = - / / e - (|fcl|+|fc2|)la ^fc 2 -fci)C(fci, fc 2 )dMfc 2 . 
-OO ^ — OO •/ — OO 

(2.45) 

Since 

r+oo 

/ dye-Vy^e' 2 ^ = T(2s)(l + 2|fc|)- 2s , 
Jo 

we conclude that 

/*°° _ /"OO 

dye-yy 28 - 1 V[E(3m(X -x-iy^^dx = -T(2s) / C(fc, fc)(l+2|fc|)- 2s dfc. 



2 

(2.46) 

The bound on V[E(/(X)|J r )] in Proposition O now follows from (|2~l3l and (|2~151 
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To make the steps (12. 45112. 46]) rigorous, we first restrict integration in (|2.45p to 
[—A, A], and then let A ->• oo. It follows from (|2.44|) that 
f A 

/ V[E(3m(X - x - iy)- 1 ^)^ (2.47) 

J -A 



-A 

POO poo 



= I H r e -(M+M)lal sin (f (** - fc i)) C{kli k 2 )dhdk 2 . (2.48) 

2 J^oo J.oo to - «l 

Multiplying (|2.48p by e _y y 2s_1 and integrating over y g (0, +oo), we obtain 

/>+oo />A 

/ dye-Vy 23 - 1 dxV[E(3m{X - x - iy)- 1 \ F)] (2.49) 

JO J -A 

= \v{2s) r (1 + \h\ + \k 2 \)- 2s Sm ^ k2 ~ h)) C(k 1 ,k 2 )dk 1 dk 2 . (2.50) 
2 J^oo J-oo k 2 - to 

We note that the integrand in (12.501) is absolutely integrable over R 2 for s > h, so 
the last step is justified by the Foubini theorem. Since E\X\ < oo, it follows from 
(|2.38l) that the kernel C(to, to) has bounded continuous first partial derivatives (see 
Lemma |2~3|) below. We split the integral in (|2.50j) into two, over 

S :={(*!, fca) : \k 2 -h\ <A-'} 

and over R 2 \ S. For (to, k 2 ) 6 S, we use 

|C(to,to)(l + |to| + |toir 2s -C(to,to)(l + 2|toir 2s | < const\k 2 -k 1 \(l + 2\k 1 \)- 2s 
which implies that the integral over S equals to 

ir(2s) / dt Sm ^ [ dkC(k,k)(l + 2\k\)- 2s + o(l). (2.51) 



2 t 



where we made the change of variables (fci, to) — > (t = k 2 — ki, k = ki). 

To estimate the integral over R 2 \ S, we restrict our attention to the quadrant 
to > 0, to > 0. The other three cases are similar. Denote Ci(t,u) = C(to,to), 
where it = to + to and t — k 2 ~ k\. We have to estimate 



du / ^(l + u)- 2s Sm( ^ Ci(t,u). (2.52) 
Integrating by parts with respect to t, we obtain 

, C 1 (t,u)dt= — *Y 7- g i(*»»)^ ( 2 - 53 ) 



cos(Ai) ,„ 

^Ci(*.«)L-.- ( 2 - 54 ) 

It is not difficult to see that the r.h.s. of (|2.53|) is bounded in absolute value 
by const^j (| log u\ + A e ) and (j2.54[) is bounded in absolute value by const A~ 1+e . 
Therefore, ther integral over R 2 \ S goes to zero as A — > oo. 
Finally, we note that the term in (|2.51|) converges to 

/oo 
C{k 1 k){l + 2\k\y 2s dk. (2.55) 
-oo 

This finishes the proof of Proposition 12.21 modulo Lemma T2.3I below. □ 
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In the proof of Proposition l2~2"l we used the fact that C(k\, fe), denned in (|2.38[) . 
has continuous bounded first partial derivatives. This is the statement of the fol- 
lowing lemma. 

Lemma 2.3. Le£E|X| < oo, andC{ki,k^) be defined as in 12.38)) . Then C(k\,k2) 
has continuous bounded first partial derivatives. 

Proof. We recall that 

C(h,k 2 ) = Cov (E(e lklX \F),E(e lk2X \F)) 

= E (E{e lklX \Fy&{e lk * x \F)) - E{e lklX )E{e lklX ). 

It follows from the Lebesgue dominated convergence theorem (for conditional ex- 
pectations), that 

dC ^k 2 ) = . Cm (mxe ik i x \T)Me ik * x m) . (2.56) 
ok\ 

Applying the Lebesgue dominated convergence theorem one more time, we obtain 
that dc ^iM) j bounded continuous function. □ 

aki 

3. Mathematical Expectation and Variance of Resolvent Entries 

This section is devoted to the estimates of the mathematical expectation and 
the variance of the resolvent entries. For z £ C \ R, we denote the resolvent of Xm 

by 

R N {z) :={zI N -X N )-\ (3.1) 

If it does not lead to ambiguity, we will use the shorthand notation Rij(z) for 
(R,N(z))ij, 1 < i,j < N. We start with the following proposition. 

Proposition 3.1. Let Xn = TT^Ww a random real symmetric (Hermitian) 
Wigner matrix \l-4\ ) (respectively il.SWl.lty) ). Then 

1 

\3mz\ e N 

uniformly on bounded subsets of C \ R, 

ERij(z) = O ^(l 3 ^ 2 ! — V\ ) l<i^j<N, uniformly on C \ R, (3.3) 
VRijiz) = O ^ML2^£J — 1^ ( l<i,j<N, uniformly on C \ R. (3.4) 

where we denote by Pi(x), I > 1, a polynomial of degree I with fixed positive coeffi- 
cients. 

If in addition, 

sup ^\(W N hf < oo, su P E\(W N ) u \ 3 < oo, 



R ii (z) = 9a (z) + 0[ ,„_, 6 . r ), (3.2) 



th 



en 



E/?:,( : ) = ()[ " )• !': ' - ; ./ ' A', uniformly on € \ ?.. (3.r>) 
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This proposition is the extension of Proposition 3.1 in 23 to the non-i.i.d. case. 
Since the proofs of (|3.2II3.5|) are very similar to the proofs given in Proposition 3.1 
in Section 2 of [23 , we leave the details to the reader. 

The next proposition is instrumental in extending Theorem ll.3l to the test func- 
tions from H s for s > 3. Our goal is to obtain an upper bound on V[(i?jv)jj(z)] 
which is integrable with respect to x — Uie z over the real line for 3m z ^ 0. 

Proposition 3.2. Let Xn — ~^^n be a random real symmetric Wigner matrix 
il.l]) , \1.4-\ ) such that the condition \1.3b}) is satisfied for some fixed m > 1 . Then 
there exists a random real symmetric Wigner matrix Tjy and a non-random positive 
sequence ejy — > as A — > oo such that the properties \2.21$2~2§^i from Lemma \2.2\ 
are satisfied and, in addition, 



IS 



(3.6) 

1 < i < m, 1 < j < N, uniformly on C \ M, where Gn(z) := (^zIn — ;^=Tjv^ 
the resolvent of — =TV. 

An equivalent results holds in the Hermitian case. 

Proof. The existence of a Wigner random matrix TV that satisfies (|2. 24112. 21?]) 
follows from Lemma 12.21 All is left to us is to show that (|3.6[) holds. Since 
P(Xjv = Tjv) — > 1 as N —> oo, we can assume, without loss of generality, that 
Tn = A^v- 

Let I be a positive constant that will be later chosen to be sufficiently large 
depending on a, a\, and 7714. We note that if 

' > T (3-7) 



\3mz\ 4 N ~ L 
then 

V[(Jto )«(*)] < n\RN(z)\\ 2 < ^ff 2 - (3-8) 

Thus, (0771) implies (|3T6|) . 
Now, let us assume that 

11 

|am*| 4 A < I" (3J) 

One can rewrite (13.91) as 



3m2 l > jyT/4- ( 3A °) 



L l/4 

1 3mz\ > 

Let us fix 1 < i,j < m. Then 

zERij(z) = 5ij + ^E(X ik R kj (z)). (3.11 



AT 



fc=l 



To estimate M(XikRkj(z)), we use the decoupling formula (see e.g. (i) in Section 
2 in [TB] and Proposition 3.1 in [T9~]). Let £ be a real random variable with p + 2 
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finite moments, and <j> a real- valued function with p + 1 continuous and bounded 
derivatives. Then 

EmO) = ]T ^lE(<^)(0) + e p+1 , (3.12) 

where the cumulants of £, 

|e p+1 | <Csnp\^+ 1 \t)\E(\^+ 2 ), (3.13) 
t 

and C depends only on p. Moreover, as follows from the proof of Proposition 3.1 in 
[19j . if supp(£) C [— K, K] then the supremum on the r.h.s. of (|3.13|) can be taken 
over te[-K,K}. 

The derivative of Rki with respect to X pq , for p ^ q is given by 

dR 

7 7T^ L = RkpRql + RkqRpl- (3.14) 
OAp q 



For p — q the derivative is given by 

dRki 
dXpp 



= RkpR P i- (3.15) 



Applying (|3. 12113.15]) to the term E(XikRkj) in (j3.11[) . we obtain the following 
Master equation 

2 



zER tJ (z) = Sa + (J 2 E[R l3 {z)tr N R N {z)] + — E[(i? w (z) 2 y (3.16) 

2(T 2 

- — E[R u (z)R t:j (z)} + r N (3.17) 



- % + a 2 E[R t] (z)tT N R N (z)} +r N + ( 3^E^)L \ , ( 3 .1 8 ) 

where contains the third cumulant term corresponding to p = 2 in (|3.12p . and 
the error due to the truncation of the decoupling formula (|3.12[) at p = 2. For k = i, 
we truncate the decoupling formula (|3.12p at p = 0. 
We will need the following lemma. 

Lemma 3.1. The following two bounds hold. 

c^j^w, *»*,(,» - o ( MjjfflWjg ) , (3 . 19) 

uniformly in z € C \ R. 

r ,v . ( P '" J -'7 E » fi "'"" a ) , (3^20) 

uniformly in z satisfying LS.10\) . where L is an arbitrary fixed positive number. 

Proof. The bound (|3.19|) follows from the first of the two bounds on the variance of 
the trace of the resolvent in Proposition 2 of [33] . It should be mentioned that the 
bound is valid provided the second moments of the diagonal entries are uniformly 
bounded and the fourth moments of the off-diagonal entries are also uniformly 
bounded ([27]). 

To prove the bound ()3.20p . one has to study the third cumulant term that cor- 
responds to p — 2 in the decoupling formula (|3.12j) for k ^ i and the error terms 
due to the truncation of (|3.12[) at p = 2 for k ^ i and at p = for k — i. 
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The third cumulant term gives 
^^[4E(^ K 3 ({W N ) lk )R lJ R lk R kk )+2E{J2 ^{{W N ) lk )R vl R kk R kj ) 

k:k^i k:k^i 

+ 2E( K 3 ((W N ) lk )(R lk ) 2 Rj k )], 

k:k^i 

where «3((Wjv)ifc) denotes the third cumulant of (PTjv)ife- Since |«3 (( Wat ) | < 
const (1714) , 

E i^i 2 ^ n^(*)ii 2 . and i^iw < ii^wii < ( 3 - 21 ) 

one observes that the third cumulant term can be bounded in absolute value by 

( n\R N (z)\r 

V \Jmz\N 

To estimate the error term due to the truncation of (|3 . 1 2[) at p = 2 for k ^ i, we 
have to consider finitely many sums of the following form 

N~ 2 E I J2 Bupl^H^II^II^I ] , (3-22) 

\k:k^i J 

where a,b,c,d,e,f,p,q,s £ {i, k,j}, the supremum in (|3.22[> is considered over 
all possible resolvents = (z — X^)~ l , I = 1, ...4 of rank two perturbations 
X$ = X N + xE lk of X N with {E ik ) jh = 5ijS kh + S ih S kj . Since 

\X lk \ < e N N~ 1/4 , k^i, e N -> as N -> oo, 

by (|2.29p . we can restrict x in the supremum in (|3 . 22[) to |x| < e^N^ 1 / 4 . Then 

R%\z) = {zI N - X^)- 1 = {zI N ~X N + xE ik )~ l = (I N + R N {z)xE lk )- 1 R N {z). 

Since by taking into account (|3.10[) 

1 AT 1 / 4 
\\R N (z)xE ik \\ < reNN- 1 / 4 < --—eNN- 1 / 4 = ^L^l 4 = o(l), 



we have 



II^WII< 11^(^)11(1 + 0(1)), 



E[|-R. w (2)l| 



and we obtain that the expression in (|3.22l) can be bounded from above by O ( x 
It follows from 

11^)11 = a- h l iy V\ ~ 1 Jm («)l _1 - (3-23) 
dist(z, bp(Xpf)) 

that one can write the upper bound as O f ^j^p^y ) • 

To estimate the error term due to the truncation of (|3.12[) at p = for k = i, 
one proceeds in a similar manner. Lemma |3. II is proven □ 

The rest of the proof of Proposition 13.21 is similar to the proof of (3.3) in 23^. 
The details are left to the reader. □ 
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4. Proof of Theorem [TT3l 

The goal of this Section is to prove Theorem 11.31 

First, we extend the estimates of Proposition 13. II to a sufficiently wide class of 
test function by using Helffer-Sjostrand functional calculus ([II], [10]) as in |23) . 
Let / € C' +1 (R) decay at infinity sufficiently fast. Then, one can write 

m N ) = --J c g=R N (z)d xd y , _:= 5 (_ + ,_j (4.1) 

where: 

i) z = x + iy with x, y € R; 

ii) f(z) is the extension of the function / defined as follows 

n=0 

here cr e C°°(R) is a nonnegative function equal to 1 for \y\ < 1/2 and 
equal to zero for \y\ > 1. 

Using the definition of / (see (|4.2|l ) one can calculate 

I - Kl + 'g) <«» 

n=0 y 

and derive the crucial bound 

?L(x + iy)\ <Canstmax(\^(x)\, l<j<l + l) \y\ l ■ (4.5) 

Directly following the calculations in Section 3 of [23 , one obtains the following 
extention to a non-i.i.d. setting of Proposition 1.1 in [23 . 

Proposition 4.1. Let Xn — T^^v be a random real symmetric (Hermitian) 
Wigner matrix \l-4\ ) (respectively H.8W1. 1(J\) . Then the following holds. 

(i) Let L be some positive number, f G C 7 (R) with compact support, and 
supp(f) C [— L,+L]. Then there exists a constant Const(L,o,o~\,m$) such that 



2a 



1 HT~2 2 -1 \ UT A\f\\c 7 ([-L,L] 

— - v 4cr z — x z dx\ < Const{L,o-,o~i,m4) — 



\E(f(Xff)ii)- fix)- — 5- v4cr 2 - x 2 dx < Const{L,o,ai,m A ) 

J-1a 27rcr N 

(4.6) 

1 < i < N. 

(ii) Let f € C 8 (R), then there exists a constant Const(a,a\,m^) such that 

\E(f(X N ) u ) - f /(sjJLv^-a^dzl (4.7) 



2fT 



< Const(a, o-i, m 4 ) , 1 < i < N. 



where ||/||n,l,+ * s defined in 
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(Hi) Let f £ C 6 (R), then 



\E(f(X N ) jk )\ < Const^aum^Mhl, l<j<k<N, (4.8) 



where ||/|| n ,i * s defined in \1.2l 
(iv) Let f e C 4 (tt), then 



II fll 2 

W( X N)ij) < Consti^a^m^^j^, 1 < i,j < N. (4.9) 



(v)If 

and f € C 10 (R), then one can improve C731), namely 



sup E\(WN) tJ \ 5 < oo, supE|(W^k| 3 < oo, 



\E(f(X N ) jk )\ < ConstM^, l<j<k<N, (4.10) 

where Const depends on sup^-^r E|(WAr)ij| 5 , and sup i N E,\(Wn)u\ ■ 
The next proposition is a corollary of Propositions 12.21 and 13.21 



Proposition 4.2. Let Xn — ~^^^ N ^ e a ran ^ om rea ^ symmetric Wigner matrix 
1 1.1]) . {l.J$ such that il.36\) is satisfied for some fixed m > 1. Then there exists 
a random real symmetric Wigner matrix Jjv and a non-random positive sequence 
ejv — ► as N — > oo such that the properties \2.21$2~2§\l from Lemma \2.2\ are 
satisfied. In addition, for s > 3, there exists a constant const s that depends on 
s, <7, <7i , and ?/i4 swc/i t/wrf for f £ T-L s 

V[/(7V%/A%] < constj^, 1 < » < m, 1 < j < JV. (4.11) 

Proof. The existence of random real symmetric Wigner matrix Tjv satisfying (|2.24t - 
I2.29[) has been proven in Lemma T2.2I Since ¥(Xn = Tjy) — > 1 as iV — > oo, we can 
assume without loss of generality that Tjy = Xjv- 

Let us first consider the diagonal case i = j. Without loss of generality, one can 
assume i = 1. Define a random spectral measure 

N 



IM(dx,u) :=^(z-A ; )|<Ml)| 2 , 



i=i 



where A; , 1 < I < N, are the eigenvalues of Xn and (j>i , 1 < I < N, are the 
corresponding normalized eigenvectors. Since by the result by Latala |17) 

supE||Xat|| < oo, 

N 

we have 

E J \x\n(dx,u) = E(|Jfjv|)n < oo, 
one can apply Proposition 12.21 and obtain 

p OO p OO 

V[/(Xjv)ii] < Const s \\f\\ 2 s / dye'Vy 28 - 1 V[{R N (x + iy))u)dx. (4.12) 

JO J-oo 
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To estimate the integral V[(Rn(x + iy))n]dx in (|4.12[) . one uses the upper 
bound (|3.6[) in Proposition 13.21 to obtain 



4[U j E / \\R N (x + iy)fdx (4.13) 

^* J —OO 

+ ^S^2 E \\R N (x + iy)f' 2 dx. (4.14) 

We will treat the first term (|4.13l) . The second term (|4.14l) can be treated in a 
similar fashion. For x € [— ||-Xjv||, +||Xjv||], we use the trivial bound 

\\R N (x + iy)\\ 2 < \. 

y 

For |a;| > ||Xjv||, we write 

1 



\R N {x + iyW < 



(x- 11^11)2 + 2,2 



Thus, 



\R N ( x + l y)\\ 2 dx<^^ + -. ,+ !:,. 



y 2 y 



Since (pi]) 

SUpE|jXAr|| < OO, 

N 

we obtain 

V[f(X N)ll ] < Const S M1 r ^ V- '/',(,/ i; + . (4.16) 

If s > 3, the integral in (|4. 16[) converges. 

In the off-diagonal case i + 1 j, one can consider the (complex-valued) measure 

N 

H(dx,u) :=22S(x- \i)<j>i(i)<f>i(j), 
i=i 

write it as a linear combination of probability measures, and apply Proposition ^. 21 
to each probability measure in the linear combination. Proposition 14.21 is proven. 

□ 



Now, we are ready to prove Theorem 11.31 Let to be a fixed positive integer. 
Denote by the to x to upper- left corner submatrix of Wn, and by R^\z) 

the to x to upper-left corner of the resolvent matrix Rn(z). Our next step is to 
compute the limiting distribution of the normalized entries of R ™ (z) in the limit 
N — > oo. In the i.i.d. sertting, this was done in Theorem 1.1 (real symmetric case) 
and Theorem 1.5 (Hermitian case) in [23) . Below, we extend these results to the 
non-i.i.d. setting. We start with the real symmetric case. Define 

T N (z) := VN (R (m) (z) - g a (z)I m y z € C \ [-2a, 2a], (4.17) 
<M*) ■= Tn(z) - gl(z)W { N m) = (r^\z) - g a {z)I m ) - &{z)Wfr\ (4.18) 
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Clcraly, T n(z) and ^n(z) are random function with values in the space complex 
symmetric m x m matrices, (real symmetric m x m matrices for real x). Define 

f 2a 1 1 1 / f a a (w)- 9a (z) . . f , 

Lp{z,w):= — jV4cr 2 - x 2 da; = <^ ««-* ^ ' 

,/-2<7 z - ^ W - X 2tT(J z I -ffa\ z ) ■ $ W = Z. 



(4.19) 

l l 



for z,w e C \ [—2cr. 2cr]. One can write ip(z.w) = E {^j^-^^jj > where 77 is a 
Wigner semicircle law (11.121) random variable. Let 

/2ct j 11 
me ffie — ^ v 7 ^ 2 - x 2 da; (4.20) 
_ 2cr z — a: w — x 2na z 

= J w ) + <P& w) + tp(z, w) + ip(z, w)) , 

/la ^ I \ 
3m 3m — 7 \/4(t 2 - x 2 dx (4.21) 
_ 2(J z — x w — x 2ito l 

= -\ {^{ziw) +(p(z,w) - (p(z,w) - <fi(z,w)) , 

/2a ^ \ \ 
mz 3m \/ 4ct 2 - x 2 dx (4.22) 
_ 2a Z — X W — X ZTTa z 

1 

= (<p(z, w) + ip(z, w) - ip(z, w) - tp(z, w)) . 

Theorem 4.1. Let Xm — ~jx^ N ^ e a ran ^ om rea ^ symmetric Wigner matrix 
fZHP , jl-4]) - Let m be a fixed positive integer and assume that for 1 < i < m 
the conditions HI .36}) and U.38\) are satisfied. Also assume that the Lindeberg type 
condition for the fourth moments of the off-diagonal entries and the Lindeberg 
type condition HI. 5) ) for the second moments of the diagonal entries are satisfied. 

Then the random field ^n(z) in _?<$[ ) converges in finite- dimensional distribu- 
tions to the random field 

*{z) = g*(z)Y(z), (4.23) 
where Y{z) = (Yij(z)) ,Y^{z) = Yji{z), 1 < i,j < m, is the Gaussian random field 
such that 

Cov(mcY u (z),meY u (w)) = n 4 (i)me g a (z)Vlc g a {w) + 2a 4 tp ++ (z,w), (4.24) 

Cov(3mYu(z), 3m Yii(w)) — n±{i) 3mg a (z) 3mg a (w) + 2a 4 tp (z,w), (4.25) 

Cov(meY u (z),3mY u (w)) = K 4 (i) fHeg CT (z) 3mg a (w) + 2a 4 ip+-(z, w), (4.26) 

Cov{fteY ij (z),VkYi j {w))=a 4 -<p ++ (z,w), 1+3, (4.27) 

Cov(3mY ij (z),3mY ij (w)) =a i ip— (z,w), i^j, (4.28) 

Cov(mcY lj (z),3mY lj (w)) = a 4 ip+-(z,w), i ^ j, (4.29) 

where ^4(7) = 7714(1) — 3<7 4 , 1 < i < m, and 7774(1) is defined in hi. 3 8} ). 

In addition, for any finite r > 1, the entries Yi l j l (z;), 1 < %i < ji < m, 1 < I < r, 
are independent if for any 1 < ?i 7M2 < r one has (ijuJii) 7^ (ii 2 ,ji 2 )- 

Now, we consider the Hermitian case. As before, we define by (14. 18)) the matrix- 
valued random field ^n(z), z € C \ [— 2<r, 2cr]. ^E'Ar(x) is Hermitian for real x and, 
more generally, &n{z) = ^n(z)* ■ 
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Theorem 4.2. Let Xn = ~^Wn be a random real Hermitian Wigner matrix A 1. Si - 
ll. 10]) . Let m be a fixed positive integer and assume that for 1 < i < m the conditions 
hi. 3b}) and hi. 3 8}) are satisfied. Also assume that the Lindeberg type condition hl.2() 
for the fourth moments of the off-diagonal entries and the Lindeberg type condition 
(1.5\) for the second moments of the diagonal entries are satisfied. 

Then the random field ^>n(z) converges in finite- dimensional distributions to the 
random field 

*(z) = g 2 a (z)Y(z), (4.30) 
where Y{z) = (Yij(z)) , 1 < i, j < m, is the Gaussian random field such that 



Cov(UKeY u (z) 




= «4 (i) 9te g a (z) *He g a (to) + 


4 

a <^ +4 


-(z,w), 


(4.31) 


Cov(3mY u (z] 


),3mY u (w)) 


= K4,(i) 3mg a (z) 3mg a (w) 4 


4 

- a tp- 


-(z,w), 


(4.32) 


Cov(meY H (z) 


>,3mY u (w)) 


= k 4 («) <Reg a (z) 3mg a (w) + 


4 

a ip+. 


-(z,w), 


(4.33) 


Cov(mtYij(z] 




= \v A {¥++{z, w) + <p—(z, 


w)), i 




(4.34) 


Cov(3m Yij (z), 3m (w)) 


= ^ i i l P++(z,w) +<p— (z. 


,w)), i 




(4.35) 


Cov(mtYij(z] 


),3mY lJ {w)) 


= 7j° 4 ( l P+-(z,w) - (p+-(w 






(4.36) 



where K 4 (i) = 7714(1) — 2a 4 , 1 < i < m, and 7774(7) is defined in HI. 38]) . 

In addition, for any finite r > 1, the entries Yi l j l (z{), 1 < i; < ji < m, 1 < I < r, 
are independent provided {ii^j^) 7^ (ii 2) ii 2 ) for 1 < l± ^ 1% < r. 

Remark 4.1. If the distribution of the entries of Wn does not depend on N, the 
random field 

T N (z)^VN(R { - m \z)-g lJ {z)I n ) , z e C \ [-2a, 2a] 

converges in finite- dimensional distributions to g^(z)(Y(z) + W^ 71 '), where Y{z) is 
independent from W^ m \ 

Below, we sketch the proof of Theorem 14.11 The proof in the Hermitian case is 
very similar. 

Proof. As in [23], one can write 

R { ™\z) = (zl m - XW - M*RMy l = (zl m - ^=w4 m) ~ M * RM 

(4.37) 

where XjtJ 1 is the 777 x 777 upper-left corner submatrix of X N , X^ N ~ m ^ is the 
(N — m) x (N — 777) lower-right corner submatrix of Xpj, 

R N ( Z ) = (zi N - m - x^-^y 1 , 

is the resolvent of X^ N ~ m ) : and M is the the (N — m) x 777 lower-left corner 
submatrix of Xn- We will denote by . . . , € R""" 1 the (column) vectors 
that form M, and by M* the adjoint matrix of M. 

It follows from Proposition 12.11 that Rn{z) is well defined for any fixed z 6 
C \ [—2a, 2a] with probability going to 1. 
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Define the m x to matrix Tn{z) as 
(T N ) %3 (z) = {W N )ij + (^{x^,R n (z)x^) - a 2 g a (z)5 t ^ , 1 < i,j < m. (4.38) 



Y N {z) = W { ™ ] +Y N (z), (4.39) 



Then 
where 

(Yjv(z))y = Yij(z) = VN^\R(z)x {l) ) - a 2 g a (z)6 t ^ , 1 < i,j < to. (4.40) 
Equations P~3Tj) and (|435)l imply 

RW(z)=(*I m - J= T N (z)) . (4.41) 

It will follow from the Central Limit Theorem for Quadratic Forms (see discussion 
below and the Appendix) that ||rV(z)|| is bounded in probability. This would imply 
that 

T N (z) = ^N^ m \z) - g„(z)I m ) = gl(z)T N (z) + o(l), (4.42) 
in probability (meaning that the error term goes to zero in probability), and 

* N (z) = VN(R( m \z)-g a (z)I m ) ~gl(z)W { N m} =gl(z)Y N (z) + o(l), (4.43) 
in probability. 

To estimate || Tjv (^) || , where Tjf(z) — Wjy + Yp?(z), we note that for fixed 

m, llW^jv™^!! is bounded in probability. Let us consider in more detail Yn(z). Assume 
that z is fixed and 3m z 7^ 0. It follows from 

EY N {z) = VNa 2 (g N (z)-g a {z))I m , 

and Proposition 13. II that EYjv(z) — > 0. Thus, 

Y N (z)ij = VN ((x {l \ R(z)x u) ) -E(a; (4) ,^(z)x (j) )) +o(l), l<i,j<m. (4.44) 

We note that the vectors x^ , 1 < i < m, are independent from R(z). In the 
Appendix, we point out that the Central Limit Theorem for Quadratic Forms also 
holds in the non-i.i.d. case under the conditions on the entries of x^ l \ 1 < i < m, 
that are equivalent to (11.36[) . This implies that ||Yjv(z)|| is bounded in probability, 
and therefore 1 1 T jv (-z) 1 1 is bounded in probability as well, which implies (14.42114. 43l) . 

To study the finite-dimensional distributions of Yfj(z), we fix a positive in- 
teger p > 1, and consider zi,...,z p € C \ K. Taking into account (|4.44l) . the 
problem is reduced to the question about the joint distribution of the entries 
y/N[(R N (zi)) llJl - E(R N (zi)) iltjl ], 1 < n < ji < to, 1 < I < p. To this end, 
we apply Theorem IA.4I in the Appendix with r — m, and 

v 

B N = E ( a S ^M z i)) + b sl 3m(R(zi))) , 1 < s < t < m, (4.45) 
i=i 

where a ( p t , b ( P t , 1 < s < i < to, l<l < p, are arbitrary real numbers, and 



(s) VN u , , 
Vn = — x > 1 < s < TO. 
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The condition (i) of Theorem IA.4I is equivalent to (|1 .36|) . The condition (ii) is 
automatically satisfied as long as 3m zi ^ 0, 1 < I < m. Conditions (iii) and (iv) 
are equivalent to 



Tr f mz(R(z)) fRt(R(w))) -> tp++(z, w), (4.46) 

— m V / 

- — Tr hm(R(z)) 3m.(R(w))) <p__(z,w), (4.47) 

— m \ / 

- — Tr f me(R(z) 3m(R(w)) -> ip+- (z,w), (4.48) 

— m V / 

1 N 

(4.49) 

1 w 

j—m-\-l 

(4.50) 

1 * 



j=m+l 

(4.51) 

for z, w e C \ [— 2cr, 2ct], 1 < i < m, where ^++(2:, ic), 9? (z,ui), and <^ (2,10) 

are defined in (|4.20M.22|) . and the convergence is in probability. To make the 
formulas (|4.49M.5lj) look less cumbersome, we label the diagonal entries of the 
(N - m) x (N — m) matrices D\z(R(z)), 3m(R(z)) by index j = m + 1, . . . , N. 

The conditions (|4.46I4.48[) follow from the semicircle law, and (|4.49H4.5l[) fol- 
low from the estimates (|3.2p and (13.4[) in Proposition 13.11 The details are left to 
the reader. Theorem IA.4I now implies that Yn(z) converges in finite-dimensional 
distributions to Y(z) for 3m z ^ 0. For z € K \ [—2a, 2a], one can replace R(z) 
by h(X^)R(z), where h satisfies (| 1 . 1 7[) and repeat the arguments above since 
F(R(z) ^ h(X N )R(z)) ->■ as N -> 00. □ 



To complete the proof of Theorem II .31 we first restrict our attention to the 
four time continuously differentiable test functions with compact support. Let 
/ G Cp(R). It follows from Theorem 11.11 and Proposition 14.11 that the result of 
Theorem 1 1 . 31 holds for finite linear combinations 
k 

^aMx) , Zl # [-2a,2a], \<l<k, (4.52) 

zj — x 

1=1 

where hi € C£°(R), 1 < I < k, satisfies (| 1 . 17|) . By Stone- Weierstrass theorem (see 
e.g [25]), one can approximate an arbitrary C^(R) by functions of the form f|4.52|) . 
Moreover, if supp(f) c [—A, A], one can choose the approximating sequance in such 
a way that supp(hi) C [—^4 — 1,A+ 1], Applying the bound (|4.9p in Proposition 
14.11 we show that 

k 



v[Vjv(/(^)« -J2 a i( h i( x N)RNh)] 



=1 
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can be made arbitrary small uniformly in N, which finishes the proof for / G C 4 (M) . 

To extend the proofs to the case of / € H s , for some s > 3, we use the estimate 
(|4.11[) in Proposition 03] and approximate such / by a sequence {f n }n>i such that 

||/ - /„|| s -> 0, as n -> oo, /„ e C C 4 (M), n > 1. (4.53) 

This finishes the proof of Theorem 11.31 

Appendix A. Central Limit Theorem for Quadratic Forms 

The appendix is devoted to the formulation of the CLT type results for the 
quadratic forms y* N By^ where yrj is a random N- vector that contains independent 
entries with finite fourth moment and B is a random N x N Hermitian matrix. 
The formulated results and their proofs are similar to the results in [4], [9] (see 
the appendix by Baik and Silverstein), and [7] since the arguments presented there 
work with small changes in the non-i.i.d. setting as well. 

First we present the case where the entries of Y/v are complex and then the case 
where the entries are real. 

Theorem A.l (Central Limit Theorem for Quadratic Forms). Let B = (&ij)i<ij<jv 
be a N x N random Hermitian matrix and y^ — {yNj)i<j<N oe & n independent 
vector of size N which contains independent complex standardized entries such that 
sup N j E | j/ jv j | 4 = TO4 < oo and E(yjy ■) = 0. Assume that 

(i) for all e > 0, 

1 N 



as N — > oo, 

(ii) there exists a constant a > (not depending on N) such that \\B\\ < a, 

N 



(Hi) jrTrB 2 converges in probability to a number 02, 



fi v ) ~k y^-i l>ii K 4(yNi) converges in probability to a number a\, 
where 

Ki{y Nl ) := E\y Nl \ 4 -2, 1 < i < N. (A.2) 
Then the random variable -j=(y* N By^— TrB) converges in distribution to a Gauss- 
ian random variable with mean zero and variance 

v 2 — ax + a 2 . 

Theorem A.2 (Central Limit Theorem for Real Quadractic Forms). Let B = 

{bij)i<i,j<N be a N x N random real symmetric matrix and y^ = {yNj)i<j<N 
be an independent vector of size N which contains independent real standardized 
entries with sup N 3 ;E|yjyj , | 4 = m 4 < 00. Assume that conditions (i)-(iv) hold as in 
Theorem \A.1\ with 

K 4 (y m ) := E\y m \ 4 - 3, 1 < * < N. (A.3) 
Then the random variable -^{y* N ByM — TrB) converges in distribution to a 
Gaussian random variable with mean zero and variance 

v 2 = a\ + 2a2- 

Finally, we formulate the multidimensional versions of Theorems IA.1I and IA.2I 
We again consider the real and complex cases separately. 
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Theorem A. 3. Let {£> s '* : 1 < s,t < r} be a family ofNxN random matrices with 
the property that (£? s >*)* = B l ' s . Let {y$ : 1 < s < r} be a family of independent 
N -vectors with independent complex standardized entries where y^) = {y^j)i<j<N> 
swp N j E|y^-| 4 = TO4 < 00, and E[(y^) 2 ] = 0. Further assume that 
(i) for all e > 0, 



N 



N r n 2 

|y$l 2 -i 



1 * 

/v ^ 



1 {||i/Kl 2 -i|>^} 







as iV — > 00 for each 1 < s < r, 
(ii) there exists a constant a > (not depending on N) such that maxi< Sjt < r ||-B S, '|| < 
a, 

(Hi) jjTr((B s ' t )*B s ' t ) converges in probability to a number a2(s,t), 
(iv) ^J2^=i(B s ' s )ti K 4(yNi) converges in probability to a number a±(s), 
where 



Then the r x r matrix 
Gn — 



K4 (2/^)=E|y^| 4 -2, 1< 1 <N. 



t<r 



converges in distribution to an rxr Hermitian matrix G such that the linearly inde- 
pendent entries are statistically independent and 9\t(G s t), 3ta(G s t) ~ A/" (0, ^a2(s, £)) 
for s 7^ t and G ss ~ A/"(0, ai(s) + a2(s, s)). 

Theorem A. 4. Let {B s,t : 1 < s,t < r} be a family of N x N real random 
matrices with the property that (B s,t )T = B t,s . Let {y$ : 1 < s < r} be a family 

(s) 

of independent N -vectors with independent real standardized entries where y N = 
(y < Nj)i<j<N and supjyj E|y^-| 4 = TO4 < 00. Further assume that the conditions 
(i)-(iv) from Theorem \A.3\ hold with 



K4 (y^)=E|^| 4 -3, l<t<N. 



Then the rxr matrix 
Gn = 



^(fo«)- B «4 , -^™-) IS „ 1Sr 

converges in distribution to an r x r symmetric matrix G such that the linearly 
independent entries are statistically independent and G s ^ ~ A/"(0, 02(3, t)) for s =/= t 
and G SiS ~ J\f (0, ai(s) + 2a2{s, s)). 
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